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' Abstract. One major problem of phase-referencing VLBI observations 

. are phase errors due to the unknown tropospheric zenith delay at each 

' antenna. These errors degrade the quality of the phase-referenced image 

C^) . and limit the achievable astrometric accuracy. We will present and com- 
pare two independent methods to estimate the zenith delay offset at each 

r S^ ■ antenna. The zenith delay offsets can then be used to correct the data. 

Oh 1 These corrections improve the quality of the phase-referenced image and 

q ■ an astrometric accuracy of 10 //as can be achieved. With this accuracy, 

\ measurements of proper motions in the Local Group become feasible. 



1. Introduction 



Phase-referencing has become a standard technique in Very Long Baseline Inter- 
ferometry (VLBI) to measure accurate relative positions of radio sources. The 
concept of phase-referencing is based on the assumption that the phase errors of 
two sources with a small angular separation on the sky are similar. One observes 
a target source between two scans on a calibrator. Then one can interpolate the 
calibration between the two scans and apply this interpolated phase corrections 
to the phase of the target source. The phase-referenced difference phase contains 
information only about the target source structure and its position relative to 
the calibrator and contains noise from interpolation errors. 

One major problem of phase-referencing VLBI observations is the unknown 
tropospheric zenith delay at each antenna. The VLBA correlator model uses 
a seasonally averaged and latitude-dependent atmospheric model (Niell 1996), 
which can misestimate the zenith delay by a few centimeters. The excess path 
length caused by the dry troposphere to a radio signal from a source with zenith 
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angle Z is given by 

/ « 0.228^P secZ(l - 0.0013 tan 2 Z) (1) 
mb 

where Po is the total pressure at the surface (see equation 13.41 in Thomp- 
son, Moran, & Swenson 2001). This excess path length will be different for 
the calibrator and the target source since both sources usually have different 
zenith angles. Hence, a residual error will remain after phase-referencing. These 
residual errors degrade the quality of the phase-referenced image and limit the 
achievable astrometric accuracy (e.g., Beasley <fc Conway 1995). If the correlator 
model misestimates the zenith delay Iq by 3 cm (~ 0.1 ns), one would expect 
after phase-referencing delay errors of 0.2 cm for zenith angles of ~ 60° and 
source separations of « 1°. However, these errors can be corrected if the true 
zenith delay at each antenna is known. 

2. Atmosphere-Corrected Phase-Referencing 
2.1. Phase fitting 

The fringe-phase of the target source after phase-referencing is mainly a sum of 
the position offset and the effects of the zenith delay error if the source structure 
is negligible. The phase induced by the position offset shows a simple 24 hour 
cosine behavior, while the phase errors from the zenith delay error follow a more 
complex dependence on the zenith angle according to equation 1. Because of 
the different behavior of the two contributions, it is possible to separate both 
effects and to estimate the position offset as well as the zenith delay error. 

One can compute a model phase with a position offset (two parameters) 
and a zenith delay error at each antenna (one parameter for each antenna) and 
perform a least-squares fit to the phase data. This method has been applied 
successfully by Reid et al. (1999) to measure the proper motion of Sgr A*. In 
Brunthaler et al. (2002) we applied this method to phase-referencing obser- 
vations of extragalactic H2O masers in IC 10 with respect to the background 
quasar J0027+5958 and achieved an astrometric accuracy of « 10 /xas. 

Table 1. Positions offsets of IC 10 in right ascension (Aa) and dec- 
lination (AS) relative to the phase center from the model fit to the 
phase data (A), the corrected map using the phase data (B) and the 
corrected map using the geodetic-like observations (C . 





A 


B 


C 


Date 


Act [/xasj AS [jj,as\ 


Aa |juasj AS |/tasj 


Act [/iasj AS [/LtasJ 


2002/01/12 
2002/01/17 


-4 112 
-11 134 


-1 103 
-8 125 


13 101 
1 135 


Mean 


-8±4 123±11 


-5±4 114±11 


7±6 118±17 



Second epoch observations of the H2O masers in IC 10 were made in January 
2002 and two observations were made within five days to check the accuracy and 
repeatability of our results. The data were calibrated with standard techniques 
using the AIPS software package. A priori amplitude calibration was applied 
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using system temperature measurements and standard gain curves. A fringe fit 
was performed on J0027+5958 and the solutions were applied to IC 10. 

Then we modeled the phase data of IC 10 as described above. The position 
offset for the two observations are given in Table 1 (A). The positions obtained 
from the two observations are consistent and have an rms of ~ 10 //as. The 
zenith delay offsets are typically of the order of a few centimeters. The values 
of the first observation on 2002 January 12 are given in the second column of 
Table 2. 

2.2. Phase correction 
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Figure 1. Phased-referenced image of a H2O maser spot in IC 10 
without correction a), with atmosphere correction from the fit to the 
phases b), with corrections form the geodetic- like observations c) and 
after phase self-calibration d). The contour levels start at 25 mJy and 
increase by a factor of \/2 in all images. The peak-to- noise ratios are 
77, 101, 100 and 151 respectively. 

The fitted zenith delay offsets can be used to correct the uv-data itself. 
This correction was done using a new version of the AIPS task CLCOR with 
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OPCODE='ATMO'. Fig. 1 shows the phase-referenced images of a strong H2O 
maser component in IC 10 with (b) and without (a) the corrections. The noise 
in the image is reduced while the peak flux increased in the corrected image. 
The peak-to-noise ratio in the corrected image is 101 compared to 77 in the 
uncorrected image. The positions of the maser feature in the corrected maps 
for the two observations are given in Table 1 (B) and agree within 10 /xas with 
the positions derived from the fit to the phases. The positions obtained from 
the fit to the phases and from the corrected maps are not independent and the 
agreement is not surprising. However, it is a consistency check for the method. 

2.3. Geodetic-like observations 

A second, independent method to estimate the zenith delay offset is to use 
geodetic-like observations. Here, a number of bright quasars with positions 
known to better than 1 mas were observed for 45 minutes before and after 
the phase-referencing observation. The quasars were observed at different eleva- 
tions with a frequency setup of 8 IFs spread over 450 MHz. A fringe-fit on the 
quasars yielded then a multi-band delay and a rate for each baseline and source. 
These multi-band delays and rates were then fit with a model that consisted of a 
zenith delay offset and drift at all antennas as well as a clock offset and drift at 
all antennas except the reference antenna. The parameters of the best model-fit 
are given in Table 2. The differences between the observed multi-band delay 
and the model are less than 0.1 ns. 

The zenith delay offsets obtained from the fit to the phases and from the 
geodetic-like observations are in general comparable although they can be dif- 
ferent for some antennas. The differences between the two methods are not 
surprising, since the model in the geodetic-like observations includes not only a 
zenith delay offset but also a drift, a clock offset and a clock drift. 

Table 2. Fit parameters for the observation on 2002 January 12: 
Zenith delay offset (to, p ) from the fit to the phase data , zenith delay 
offset (to) and rate (tq) and clock offset (At) and rate (At) from the 
geodetic-like observations. The errors are the formal errors of the least- 
squares fit. 



Antenna 


to, p [cm] 


r [cm] 


fo [cm h 1 ] 


At[ns] 


At[ns d" 1 ] 


BR 


2.5±0.6 


1.6±0.1 


0.4±0.03 


-0.24±0.01 


-1.52±0.03 


FD 


4.4±0.3 


1.7±0.05 


0.3±0.01 


-0.08±0.01 


-0.84±0.02 


HN 


7.0±0.2 


6.4±0.1 


0.3±0.03 


0.37±0.01 


0.47±0.04 


KP 


1.0±0.3 


-0.1±0.05 


0.3±0.01 






LA 


-0.4±0.4 


-1.7±0.07 


0.3±0.02 


-0.12±0.01 


-0.17±0.02 


MK 


-3.3±0.1 


-4.0±0.06 


0.3±0.01 


-0.93±0.01 


-0.97±0.03 


NL 


9.4±0.4 


7.8±0.1 


0.2±0.02 


0.42±0.01 


0.51±0.03 


OV 


4.7±0.4 


2.2±0.07 


0.3±0.01 


0.00±0.01 


-0.62±0.02 


PT 


-0.2±0.4 


-1.8±0.07 


0.2±0.01 


0.08±0.01 


1.05±0.02 



The fit values of the zenith delay and clock offsets can also be used to correct 
the uv-dat&. A corrected image is seen in Figure 1 (c) and the peak-to-noise ratio 
is 100. The positions of the maser feature in the corrected maps are given in 
Table 1 (C). They are also consistent with the positions of the fit to the phases. 
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3. Conclusions 

We have presented a technique to correct phase-referencing observations for 
inaccurate zenith delays in the VLBA correlator model. The quality of the 
atmosphere-corrected phase-referenced images can be improved by this tech- 
nique. The peak-to-noise ratio in the corrected images increases by ^ 30 %. By 
using this technique, an astrometric accuracy of 10 fias can be achieved. 

Two independent methods were presented to estimate the zenith delay off- 
sets at each antenna. The images obtained from the correction with the zenith 
delay offsets from the fit to the phases (Fig.l (b)) and from the geodetic-like 
observations (Fig.l (c)) are very similar and the positions measured from these 
images are consistent. This is a proof of concept for our technique. 

The astrometric accuracy of this technique enables us to detect significant 
proper motions of galaxies in the Local Group out to a distance of 800 kpc within 
a few years. One needs at least three epochs of observations to make reliable 
statements about measured proper motions. Further observations are underway 
and the results will be presented elsewhere. 
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